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Abstract 

According to theoretical physics the cosmological constant (CC) is expected to be 
much larger in magnitude than other energy densities in the universe, which is in stark 
contrast to the observed Big Bang evolution. We address this old CC problem not 
by introducing an extremely fine-tuned counterterm, but in the context of modified 
gravity in the Palatini formalism. In our model the large CC term is filtered out, and 
it does not prevent a standard cosmological evolution. We discuss the filter effect in 
the epochs of radiation and matter domination as well as in the asymptotic de Sitter 
i— | future. The final expansion rate can be much lower than inferred from the large CC 

without using a fine-tuned counterterm. Finally, we show that the CC filter works also 
in the Kottler (Schwarzschild-de Sitter) metric describing a black hole environment 
Jh with a CC compatible to the future de Sitter cosmos. 

ox) 



1 Introduction 
> 



The starting point of this work is a CC or equivalently a vacuum energy density A of 
enormous magnitude. This expectation is suggested by contributions to the CC coming 
from phase transitions in the early universe, zero-point energy in quantum field theory 
or even from quantum gravity. In general, all these parts are of different magnitude and 
probably unrelated to each other. Hence, the sum A of all terms is dominated by the 
largest contribution. Since other energy sources dilute with the expansion of the universe, 
the CC will eventually take control over the cosmos. Depending on its sign the CC would 
induce in the very early universe either a Big Crunch or an eternal de Sitter phase with 
a very high Hubble rate H oc A. Obviously, the standard Big Bang evolution does not 
happen in this case. 

The simplest way to avoid this problem is the introduction of a CC counterterm A ct , 
which makes the sum | A + A c t | smaller than the currently observed critical energy dens- 
ity PcO ~ 10~ 47 GeV . For concreteness let us assume that A ~ -Mg W were related to the 
electroweak phase transition at the energy scale M ew ~ 10 2 GeV. Then the counterterm 
must be extraordinarily close to (—A) requiring an enormous amount of fine-tuning, 



A c t 
1 + -P 
A 



< 



PcO 



A 



,. , , . -55 



W~ bb . (1) 



Apart from the fine-tuning of the classical counterterm, the situation is even more involved 
when quantum corrections are included, cf. Ref. [lj for an elaborated discussion in the 



context of the electroweak sector of the standard model of particles. Moreover, the problem 
worsens when A is dominated by higher energy scales, possibly originating from grand 
unified theories where A ~ (10 16 GeV) 4 or quantum gravity with A ~ (10 19 GeV) 4 for 
instance. Summing up, the fine-tuning of the CC is considered to be one of the most severe 
problems in theoretical physics 015]. In addition, the current accelerated expansion of the 
universe |H [6] can be explained very well by a tiny CC of the same magnitude as the 
energy density of matter, giving rise to the so-called coincidence problem. For the latter 
problem many explanations have been proposed [7] [8] [9] [TO] 11] , which induce late-time 
accelerated expansion. However, most of these models tacitly assume that the large A has 
been fine-tuned away and thus they do not address the big CC problem. 

Without fine-tuning we have to accept the existence of the presumed huge CC, and 
we have to find a way to neutralise its effects in order to obtain a reasonable cosmological 
evolution. Along this line, several proposals have been made, e.g. relaxation models for a 
large CC in the context of matter with an inhomogeneous equation of state (EOS) [12], 
or in the LXCDM framework [13] with a variable cosmological term |14| 115]. see also |16| 
[T7] fTS] IT9] 120] 121] , Removing or filtering out vacuum energy has been investigated e.g. in 
Refs. ]22l|23j|2ll25], and it is a feature in unimodular gravity [26] [27] [28] [29l [30] . Recently, 
a CC relaxation model has been discussed in the context of modified gravity with an action 
functional f(R,G) involving the Ricci scalar R and the Gaufi-Bonnet invariant G in the 
metric formalism |31[ [1] , where the action is varied with respect to the metric g a b only. 

In this work, we also consider a modified gravity model with an action functional / 
in terms of the Ricci scalar R and the squared Ricci tensor Q = R a bR ab . However, here 
we apply the Palatini formalism, where the metric g a \, and the connection T^ c are treated 
independently by the variation principle. In contrast, the metric formalism requires from 
the beginning that 

KM = 29 ad (9dc,b + 9bd,c ~ gbc,d) (2) 

is the Levi-Civita connection of g a b, whereas the Palatini connection depends also on the 
matter sector. As we will show in this paper, this allows the construction of a filter for 
a large CC, and the results will be similar in effect to unimodular gravity. However, in 
our setup the CC is not eliminated completely, but it appears in suppressed corrections. 
Furthermore, we investigate the filter effect from the early universe till the asymptotic 
future in addition to black hole environments. It turns out that finite vacuum energy 
shifts originating e.g. from phase transitions can be neutralised, too. 

There is an interesting conceptual difference to the CC relaxation models in the metric 
formalism, where the large vacuum energy is not filtered out from the total energy content, 
but gravity is modified such that A does not induce large curvatures. This happens on the 
level of the Einstein equations, i.e. by solving differential equations to obtain concrete low- 
curvature solutions. In the Palatini framework of this paper, we will show that the large CC 
can be removed already in an algebraic way before solving differential equations. Moreover, 
the latter are only of second order, whereas the metric version of modified gravity generally 
involves a higher differential order, signalling the existence of new degrees of freedom, which 
can be the source of new instabilities and other problems. More differences between both 
formalisms will become visible in the forthcoming discussions in this paper. 

Some work on the Palatini formalism can be found e.g. in Refs. |32| l33| 134"] [35] 136] 137] . 
and comparisons with the metric and other formalisms were made in Refs. |38| [39] |4"0] |4"T] 
2] S3] E]. Non-trivial properties in Palatini models have been investigated in Refs. [45 1 
6] U7] US], and finally the resolution of the Big Bang singularity has been proposed in this 



context |49| [5U] , possibly in connection with loop quantum gravity |51j . 

The paper is organised as follows, in Sec. p2] we briefly reproduce how to solve f(R, Q) 
modified gravity models. In Sec. [3] we present our model which filters out the large CC. 
The results will be applied to cosmology in Sec. [4] discussing the radiation, matter and 
late-time de Sitter era. Finally, in Sec. [5] we show that the CC filter works also for the 
Kottler (Schwarzschild-de Sitter) solution describing a black hole in the presence of a CC. 
We conclude in Sec. [6] and give an outlook to future developments. 

In this work the speed of light c and the Planck constant h are set to unity, the signature 
of the metric is (— 1, +1, +1, +1). 

2 f(R, Q) modified gravity in the Palatini formalism 

The action of our setup is given by 



d x 



g\\f(R,Q) 



+ S mat [g a h(t>], (3) 



where g ab is the "physical" metric, on which the matter fields in 5 ma t propagate. The Ricci 
scalar R and the squared Ricci tensor Q depend on both the metric and the connection T bc 
while the Ricci tensor R ab is defined only in terms of the latter: 

R ab [r] = ri b>e -r e eb>a + ri b rj e -ri f r f eb (4) 

R[g,T] = R a a =g ab R ab (5) 

Q[g,T] = R ab R ab = g ac g bd R ab R cd . (6) 

Note that in general these quantities are different from their metric versions, and the 
absence of torsion implies that the connection is symmetric. Moreover, we restrict our 
discussion to the case of a symmetric Ricci tensor R ab in the action (pi). This property is 
not automatic even for symmetric connections, which was shown recently in Ref. |52] , 

With these preliminaries the variation 2 5S/5g ab = of the action functional S with 
respect to g ab yields the modified Einstein equations 

fR R m + 2 fQ R m R a ~ ^m f = T m ') ( 7 ) 

where the energy-momentum tensor T m n emerges from the term 5 ma t, and fn and /q are 
partial derivatives of / with respect to the scalars R and Q. Moreover, we obtain from 
5S/5T bc = the equation of motion (EOM) for the Palatini connection, 



V a 



g\(f R g mn + 2f Q R mn ) =0, (8) 



where V a denotes the covariant derivative in terms of the yet unknown connection T bc . 

At this point we should remark that in the metric formalism the term Q yields EOMs 
with higher-order derivatives and problematic extra degrees of freedom. Generally, it is 
difficult to avoid instabilities, e.g. of the Ostrogradski-typdj [53j . In contrast, the Palatini 

In the metric approach this type of instability can be avoided in gravity actions depending only on 
the Ricci scalar R and the Gaufi-Bonnet term G. It would be interesting to compare both approaches in 
the context of f(R, G) models. However, G contains the squared Riemann tensor, and to our knowledge 
no method to solve the corresponding Palatini EOMs has been found yet. 



formalism provides second-order EOMs for our scenario just as standard general relativity, 
and problems from extra degrees of freedom do not occur. 

In the following, we work along the lines of Ref. [39], where a procedure for solving the 
Palatini EOMs is discussed. The strategy is as follows: the formal solutions to both EOMs 
in (IT]) and (pi) relate in an algebraic way the geometrical scalars R and Q with the physical 
metric g a \, and the energy-momentum tensor T m n , respectively. Thus, one can express R 
and Q in terms of the matter content alone by solving this set of equations. Subsequently, 
the results are plugged back into the formal solutions to obtain the explicit form of the 
Palatini connection T and all the quantities derived from it. Since the results will involve 
the metric g a b and its derivatives, it is possible to relate the cosmic expansion rate with the 
matter energy density, for instance. In the rest of this section, we explain the procedure 
for general f(R, Q) models. 

First, let us write Eq. (|7j) in matrix form by introducing the 4 x 4-matrices P = R^ = 
R ma g an and T = T^, whose entries are just the components of the corresponding tensors 
components. With the identity matrix / we obtain 

f R P + 2f Q (P) 2 -^fi = f, (9) 

and the trace of this equation reads 

f R R + 2Qf Q -2f = T, (10) 

where R = tr(P), Q = tr(P 2 ) and T = T™. 

For determining the Palatini connection T we introduce the auxiliary metric h mn and 
consider the equation 

~~ 0, (11) 



v a [r] [V\h\h™ 

where V a [r] is the covariant derivative in terms of T. Consequently, in order to solve this 
equation the connection T must be compatible with h mn , i.e. it has to be the Levi-Civita 
connection of h mn , 

nM = \h ad {h dc<h + Vc - hbcd), (12) 



just as in general relativity. Hence, we convert Eq. (jSJ) into (11) by defining h mn in the 
following way, 

^\h\h- 1 = ^/\g\9- 1 ^ with t = (fj + 2f Q P) , (13) 



where the metric h mn (and analogously g mri ) has been written in matrix notation as h = 
hmn with its inverse /i -1 = h mn . Calculating the determinant of both sides, we find 
h 2 ^ 1 = h = gdet S, which allows to eliminate y/\h\ and finally yields 



h- 1 ^ 9 S - , h mn = h=^\det£\£~ 1 g. (14) 

|det t\ 



by T% c [h] in (12) if h mn is defined as in (14). 



Since the connection in (12) solves Eq. (11), the formal solution of Eq. <\8h is also given 



The remaining step is to find P which requires an explicit form for the energy-momentum 
tensor T m ™. Here, we assume that the matter sector can be described by a perfect fluid, 

T m n = (p + P)u m u n + (p - A)C, (15) 



where p and p are the energy density and pressure of (ordinary) matter, including e.g. 
dust and incoherent radiation. And u m denotes the corresponding 4-velocity vector of the 
matter field. A represents the energy density corresponding to the cosmological constant, 
and it contains all vacuum energy contributions. Thus we require p ^ —p without loss of 
generality. Next, let us write the matrix expression Mm in the following way 



(2/q) 2 M 2 
M 



x I + p u m u n 
2f Q (p-A) + f Q f + ±f R 



v ■= Vq(p + p)- 

By explicit calculation one can check that 

c a ■ 2f Q M = xf + y u m u r ' 
with 

y ■-- 



-x + c b - s/x 2 + p(u m u m ) 



(u m u r ' 



(16) 
(17) 

(18) 
(19) 

(20) 
(21) 



is a solution to Eq. (16), which yields P. The constants c a ^ = ±1 and the sign convention 2 
for x = y x 2 will be fixed later by consistency considerations. 

Now, we have to determine the scalars R and Q, which follow from the trace equa- 



tion (10) and the trace of (20) 



c a (2f Q R + 2f R ) = 4x-y, 



(22) 



where u m u m = — 1 will be used from here on. In the last equation we eliminate all roots 
by squaring twice, which results to 

(j2f Q R + 2f R f + 8x 2 + M ) 2 = 36x 2 (2f Q R + 2f R f , (23) 

where c a ^ = ±1 and odd powers of x have dropped out. Solving this algebraic equation 



together with ( 10 ) is the tough part of the Palatini formalism. Once this task is achieved, 



R, Q and P = R^ in ( 20 ) are given as functions of p, p, A only, and can be used to calculate 



the connection T. For this purpose, let us write the matrix S in (j 1 3|) as 

t = 

with L\ 



2f Q M + ^f R I = L X I + L 2 C^ = £i (i + -£■ w r 
1 



c a x + -f R , L 2 := c a y. 



(24) 
(25) 



Using det(I + a m b Tl 
reads 



1 + b m a m we find det(S) = L\(L\ — L 2 ), and the inverse matrix 



-.-l 



1 



L 



I 



Lo 



U-jjiU . 



Li{L x -L 2 ) 

2 Changing the sign of x in Eqs. (20 1 and (21 1 by x — > — x corresponds to c a 



-c a , and therefore it 



does not yield a new solution. Here, we use the convention y/x 2 + /i(u m u m ) — x^/l + /j,(u m u m ) / x 2 , and 
the second possibility — x\/l + ■ ■ ■ would just mean Cb — > —Cb- 



Finally, we have all ingredients to write down the explicit form of the auxiliary metriqj 
from Eq. (14) 



where 



van 

ran 



nig. 



Lo 



Li-L 



- ! il rn %L n 



n -i j g mn + <m u n | ^ 



J2:= 



|detE| _ Vl^(£i-L 2 ) 



L, 



U 



(26) 
(27) 

(28) 



Once h mn is known, the connection r? follows directly from Eq. (12), and subsequently 
the Ricci tensor Q and the scalars in Eqs. (pi) and (pi) can be calculated. 



3 Relaxing the CC with a filter 

In this section we study a modified gravity model to relax the CC in the Palatini f(R, Q) 
framework. Motivated by earlier work in the metric formalism |31[ [1] we consider the 
following ansatz for the gravity action 



f(R, Q) = kR + z with z:=/3 



R n 
B~™ 



B = R 2 -Q, 



(29) 



where k and ft are non-zero constant parameters and n and m positive numbers. In the 
metric formalism model [31] the structure of the function B was enforcing the universe 
to expand like a matter dominated cosmos even when the matter energy density p m was 
much smaller in magnitude than the vacuum energy density A. However, in the Palatini 
framework the function B is not known in the beginning, and it is necessary to investigate 
under which circumstances a relaxed cosmological expansion behaviour can be obtained. 
Moreover, we will see in the following that z is not a correction to the Einstein-Hilbert 
term but a crucial part of the action functional /. Therefore, one should refrain from 
considering the limit z — > 0. 



From Eq. ( 29 ) we find 



Ir 



kR + tiz z 
^ 2/ Qj R, f Q = m-, 



and the trace of the stress tensor (15) reads 



-4A + 3p - p. 



As a result, Eq. (10) provides the first equation for finding R and B (or Q), 
jz = kR — 4A + 2>p — p, 7 := (n — 2 — 2m), 



while the second one is given in Eq. (23), explicitly 

= f Q R*S- i + f 2 Q S 2 + f Q R- 2 S 1 



(30) 



(31) 



(32) 



(33) 



The relation in Eq. (261 is called a disformal transformation |54j . which is used e.g. in MOND theor- 



ies |55| and scalar field models for dark energy |56| . 



72 



kR + p + p + -z(2 + 2n + 7) 



S 2 := 4[z 2 (-17n 2 + 4n(2 + 7 ) + 4(2 + 7 ) 2 ) 

+ 12z(p + p)(2-n + 7 ) + 6zKi?(4-5n + 27) 

+ 90 + p) 2 - 9(kR) 2 ] 

Sx := 2A[z{KR + nz) 2 {n-2--f)] . 



(34) 



(35) 
(36) 



Apparently, this set of equations is quite complicated, thus we will solve it approximately. 
First, we consider in this work only the epochs when the large CC A dominates over all 
other energy sources. Then Eq. (32) implies the relation z = (3R n /B m = 0(A) suggesting 
that B/R 2 and Jq = B/(mz) are relatively small quantities, which may be used as expan- 



sion parameters in Eq. (33). At this point we cannot prove this suggestion quantitatively 



because R and B are not known yet. However, it will be confirmed later by Eqs. (43) 



and ( 45 ) . As a result of assuming that Jq is sufficiently small, the term 53 proportional 
to fn in (33) will be the most important one, and a good zero-order solution can be found 
by neglecting the other terms /q S2 + JqR~ 2 S± and solving only S3 = 0. 

Our goal is a relaxed universe, i.e. one which is not dominated by the large CC term, 
and it can be realised by requiring that the Ricci scalar R following from S3 = is 
free from large 0{z) contributions. This happens when the parameter n is restricted by 
2 + 2n + 7 = 3n — 2m = 0, which eliminates the 0{z) term in Eq. (34). We will apply this 



condition from now on, hence Eq. (33) can be written as 







kR + r 

2 ( B 
-mz —77 
9 \R 2 



(37) 



+ 



27™ \W) 



1 + (3kR + 2r) + -7 ^ 

mz 4:{mz) z 

2 / o n \ 2 



((«#)' 



r 2 ) 



1 + 



3kR 

2(mz) 



with r := p + p. 



From the first line in the last equation one clearly observes that CC terms with EOS p = —p 
do not contribute to the Ricci scalar at leading order. In other words, the CC is filtered 
out from r = p + p, which describes the matter sector. Note that z = C(A) still appears 
in suppressed correction terms. 



In the following, we will solve Eq. (37) in situations relevant for cosmology. A first order 
correction to S3 = can be found by keeping only the term z(B/R 2 ) from the second line 



in (37), which leads to 



kR + r = -mz 
9 



R 2 ) 



9 



-mz 



a 



R 



(38) 



where we used z = f3R n /B m from Eq. (29) with 3ra 
a 7 -order polynomial equation in R, 



2m in the last step. Thus, we obtain 



(kR + r) 3 (kR) 4 = L 3 



9 



-mz 



a 



(39) 



which clearly shows that R is a function of p and p only. Finally, we find the approximate 
solution 



kR 



-r + D + O 



D< 



with D 



(40) 



when the matter-related quantity r lies in the range \D\ -C r <C |A|, which we refer to as 
the early-time limit or the limit of large energy density in r. Moreover, kR < because 
r > for ordinary matter, and kR will remain negative for decreasing p because kR — > 



is not a solution of Eq. (39). 



In the opposite limit, 



0, which is called the late-time limit in the following, we 



obtain from Eq. (39) 



kR = p e - -r + O 



P< 



with p e 



T 7 



(41) 



indicating that kR approaches the negative constant p e for vanishing matter. For a given 
value of A or z the parameter j3 must be chosen adequately for L^ = p\ < 0. In Sec. 4.2 we 
will show that (—p e ) is close to the critical energy density in the asymptotic future, which 



corresponds to the tiny observed value of the effective CC in Eq. (83). The parameter 



dependence of the solution set of the full equation (37) might further constrain (3, however 
this has to be determined numerically. In Fig. [I] we show an example for m = 3 and p e < 0, 
which nicely demonstrates the validity of our approximations. Note that not all values of m 
might allow physically reasonable solutions. We will discuss some examples for f3 at the end 



of Sec. 4.2 However, in this work we concentrate on analytical results, and the complete 
parameter dependence as well as the case L3 > will be investigated elsewhere. 

Via z = (3(R 2 / 3 /B) m from Eq. (29) it is straightforward to obtain B and B/R 2 from the 
approximate solution of R. In the following we denote subdominant corrections by e <C 1. 



Accordingly, at early times Eq. ( 40 ) yields 



B 

B 
R 2 




2D 
3 r 



+ 0(e 2 ) 



(42) 
(43) 



whereas from Eq. (41) we obtain the corresponding late-time results 



B 

B 
R 2 



Pi \ g m 



-mz 



-1 



-1 



\-2 



4 r 
7 p e 



2 r 
7 p e 

+ 0{e 2 



+ 0{e 2 ) 



(44) 
(45) 



These results confirm that B/R 2 is sufficiently small to justify that we neglected some 
terms in Eq. ( |37[ ), hence, we have found consistent solutions for R and B. More support 
for the validity of the approximations is provided by the numerical solution of the complete 
Eq. (37) in Fig. [I] 

With R and B as functions of only the matter-related term r, it is possible to calcu- 
we have to define the square root of x 2 in (20). Since we work 



First, 



late L\ y 2 in (25). 

in the limit \Rr/B\ S> 1, which implies \JqR\ 

the dominant term in Eqs. 



and we choose the convention x 
Eq. Jl8| divided by {f Q R) 2 



\mzR/B\ 3> z/R, we identify /qR to be 



(30} and ( |18[). Thu s f R w 

fQRy/l H , where 

Next, we plug x « -|/r p 



-2f Q R and x 2 « \f 2 R « (f Q R) 2 , 
■ ■ denotes the remaining terms in 
fnR into the trace equation (22) 



c a {2f Q R + 2/r) = 3x + c b -xJl- (2f Q r)x 



(46) 



7 
6 - 
5 - 
KR/p e 4 
3 
2 



exact result 
early time 
late time 



012345678 

-r/p e 

Figure 1: Numerical solutions for kR < as a function of r = p with p = 0. The plot 
shows the exact result in Eq. (|37[) (black solid curve) and the approximations at early times 



(red dashed) from Eq. (40) as well as at late times (blue dashed-dotted) from Eq. (41), 
respectively. In this example with m = 3 the parameter /3 < has been chosen such that 
— |A| ~ 358 p e only for numerical reasons. A more realistic ratio of \A/p e \ would be much 
larger, however the results would not differ qualitatively, which is true also for other values 
of the EOS p/p as long as r <C |A|. In any case, kR approaches — r (dotted diagonal line) 
in the region — p e <C r <C |A|. 



which gives in leading order c a (—2)fqR = fqR{2> + q,). Obviously, c a = q, = — 1 is the 
only solution 4 ] in the large |/q-R| limit, and we find 



Li = -x+-f R 



-2f Q R + R- 1 ( mz + \kR - i r ) + |^r (1 + e) , (47) 



L 2 = -x-xJl- (2fQr)x 2 

i /4 1 \ (mz) 2 

-2f Q R + R- 1 i^-mz + -(kR + r)\ + fj-^ (1 + e 

where we used a series expansion of the root in 

2 



(48) 



f R -2f Q R + R-"[KR+-mz), (49) 



, „ / 8mz + 3KR-3r (2mz + 3kR) 2 . . 

X = fQR f W^ + 36^ • (50) 

Above and from here on e denotes small terms of the order kR/z, r/z or B / R? . Moreover, 
A has been eliminated in favour of 

(kR - 4A + 3p - p) . (51) 



6 + Am 
Finally, we write down the series expansion of 

1 mz / IB, A . 

L '- i2 = -3lf( 1 + 2i? (1 + e) )' <K) 

which appears in the auxiliary metric h a \,. Now, we have all ingredients available for 
discussing solutions to the Palatini field equations in the next sections. 

4 Cosmology 

For investigating the cosmological evolution the physical metric is of the spatially flat 
Friedmann-Lemaitre- Robertson- Walker (FLRW) type with goo = — 1 and ga = a?(t) in 
Cartesian coordinates, where ait) is the scale factor as a function of cosmological time t. 
The matter component is at rest (ti m = — 5™) in these coordinates. Accordingly, the 
auxiliary metric h a b in Eq. (|26|) is completely defined by the diagonal elements hu = Sign 



and /loo = ^(#00 — -^2/(^1 — -^2))) which follow explicitly from Eqs. (47), (48) and (52), 



y /\L x (L 1 -L 2 )\_, T ,,/X , 2D _ Wl ; ^ 



O = Vl ^ ^=sgn(L 1 )^|-(mz)2i3-i|(l+ e 2 ), (53) 

^00 = ^y zZ V = ^f-6^Vl-|,(l + £ )V (54) 



Li-L 2 V B J \ R 2 



Note that our scenario with cj = — 1 relies from the very beginning on the existence of matter, because 



otherwise the vector u m in Eq. (16 1 were absent and Cb — +1 would be required in Eq. (21 1. However 



matter is a fact of reality, which singles out q, = — 1 in our model. Consequently, Eq. Ell is the correct 
solution in the r — > limit. 

10 



Note that if B and L\ are positive, the metric h a \, has the same signature as g a i,. Moreover, 
in the following z will be treated as a time-independent constant proportional to the CC A 
since the corrections can be subsumed in the e ~ 0(kR/z) terms. At this point we are 
ready to apply the results for B and R 2 /B that we found in Eqs. (42), (43), (44) and (45), 
respectively. 



4.1 Early universe 



We begin with the epoch, where the matter energy density p (including dust and radiation) 
is well below the large CC A ~ z in magnitude but above the asymptotic future energy 
density \p e \- Therefore, our discussion will be valid for most parts of the radiation and 



matter eras. According to (53) and (54) we find 

i 

D 



n 



h 



00 



Cl 



c\ 



D 



1 + eei— + 0{e 2 ] 
6r 



c\ = const. 



4 \ -I 

-mz/ej L 3 4 D 



D 



l-ee 2 - + 0(e 2 ) 
r 



(55) 
(56) 



where r = p + p and D = (L^/r^y 1 '^' was introduced earlier in Eq. (40). The sign under 
the root in the constant c\ = \/±3(g7n£) 3 (— n) 2 may be chosen such that consistency with 



Eq. (53) is obtained. However, since c\ drops out from the connection T(h) in Eq. (12) 



we skip this question for the moment. In addition, we have introduced above the purely 
technical parameter e = 1 which counts the powers of the small quantity \D\ <C r, whereas 
ei : 2 = 1 just denote first-order correction terms. 



According to Eq. (12) the non-zero components of the (symmetric) connection read 
(withi= 1,2,3) 



1 oo 



r° 



1 i0 



\r 7 r D _ , 9 

2r br r 



a „ / a 
e D 6- 

8mz V ci 



a 



+ C(e 2 ), 
7 rD 



eei- + 0(e 2 

18 r r 



(57) 
(58) 
(59) 



For a simpler presentation we assume from now on a constant matter EOS ui = p/p > 
— 1, which yields a simple expansion law for our matter component via its conservation 
equation^ 



p + 3-p(l + uj) = 44> p cc a 



-3(l+o;) 



R 2 



Hence, from Eq. (|4j) we obtain for the non-zero components Rq 

i? 3 3 of the diagonal Ricci tensor R^ = g cb R a c the following results 



R ° and R { 



(60) 



Ro 



9 -Koo 



-(w + 3) 2(w + l) 



+ 



7 D 

~^Pm 



(ei(19w + 21) + e 2 3(u + 3)) + e 2 2- 



o 



+ 0(e 2 



(61) 



Since the matter action 5 ma t in (|3| does not involve the Palatini connection, the covariant derivat- 
ive V„ = Vn[ff] in the matter conservation equation V„T m " = contains only the Christoffel symbols of 
the physical metric g a b- 
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R; 



9 n Rn 



omz 



a \ /r, r\ a 

- (3w + 5) + - 
a / a 



+ 0(e 2 ). 



(62) 



Remember that the Palatini connection T^ c given above and all derived quantities like R 
or R a b are different from the metric versions, and they should not be compared with them. 
Observables are related to the scale factor a in the physical metric g a \,- Due to Ri = 0(e) 
the leading component of the Ricci scalar R = Rq + 3-Rj is just the 0(e°)-term in Rq: 



R=-(uj + 3) I 2(w + l) 



+ 



0(e). 



(63) 



Comparing this result with kR = —p(l + u) + 0(e) from Eq. (40), we quickly deduce the 
leading scale factor behaviour. The situation is similar to general relativity because the 
power-law ansatz a(t) oc t s , s = const, for the scale factor represents a reasonable solution. 
Equipped with the corresponding Hubble rate H = s/t and a/ a = s(s — l)/t 2 , Eq. (63) 
yields R oc t~ 2 at zero order, and p oc a~ 3 ( l+UJ > oc t^ 3s ( 1 + w ) -vvill be proportional to R if 



3(w + 1) 



in a(t) oc t s . 



Hence, we find from Eq. (63) 

2 a „ 2 1 + 3^ 



il 



3(w + l)f 



a 
a 



H 



R 



-(oj + 3) 2 H 2 + 0(e) 



implying the modified Friedmann equation 

T 2 4 w+1 



(-«)^ 



3(w + 3) 2r ' 



(64) 



"VCl + w), (65) 



(66) 



which is the main result of this section. Accordingly, a dust dominated universe with 
EOS uj = leads to 

(-k)# 2 = 4p, (dust) (67) 
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which is not very different from the radiation dominated cosmos with u) = i. In the latter 
case, the modified Friedmann equation reads 



-k)H z 



25 



p, (radiation) 



(68) 



where p is the radiation energy density. Obviously, dust matter and radiation influence 
the cosmic expansion in almost the same way as in standard general relativity if we choose 
the parameter k = — A • 3/(87rG , 7v) with Newton's constant Gn- However, comparing 
Eq. (67) with (68) indicates an increased expansion rate in the radiation era, which can 
be expressed by a higher effective Newton constant G Tax ± = 35 Gw = 1,08 Gat. This 8% 
difference seems to be well within current bounds on the variation of Gn- For instance, 
in the context of Big Bang nucleosynthesis the bound G ra d/GN = 1, 10 ± 0, 07 was given 
recently in Ref. [57J, which can be related also to constraints on (additional) relativistic 
degrees of freedom |58[ [59] I60| . Accordingly, the small difference in the expansion rates 
above is not in conflict with recent observations. 

As a more technical point we remark that k = +1/(8^(7^) is positive in general re- 



lativity, however, in our setup the negative term kR < is only one part in the action (29) 



12 



which is completed by the crucial term z. Consequently, k < does not correspond to 



a negative Newton constant in general relativity. Another interesting feature of Eq. ( 66 ) 
is that even if we had allowed additional vacuum energy (uj = — 1) contributions in p, 
it would not influence the expansion rate H because of the vanishing right-hand side of 



Eq. (66). An example for this are finite shifts \5A\ <C |A| in the vacuum energy density 



emerging from phase transitions. 

Comparing the Palatini setup with the CC relaxation models in the LXCDM frame- 
work |14l [T5] and in metric f(R, G) |31[ Q] modified gravity, we find another advantage. 
Here, it is not necessary to include a parameter which controls the transition from radiation 
to (dust) matter domination. It simply happens when the energy density pdust oc a~ 3 of 
dust matter overtakes the radiation density Eradiation oc a -4 as a result of the growing scale 
factor a. Therefore, the relation of the energy content with the expansion of the universe 
is very similar to general relativity with a negligible CC. 

For completeness we list all results including a first order correction in the scale factor 



a(t) 



t \ 3(^+1) / Y) 



(69) 



with the constant e a to be determined below. The components of the diagonal Ricci tensor 
read 



Ro 

Ri, 



(u + 3) 2 

3t 2 {l+Uj) 2 

D{uo + 3) 2 



1 + e- 



D 

P 



7[4ei(5 + 4u) + 3(3 + w)(e 2 + e a (23 + 19w))] 



3(3 + u ) 



+ G(e 2 



+ 0(e 2 ), 



Amzt 2 (l + u) 2 
yielding the Ricci scalar R = Rq + 3R{ and the squared Ricci tensor Q, 



R 

Q 



Ro + 



3t 2 (l + uj) 2 \ + 4mz { ' 



R a b R b a = (RoY + 3{R t Y = (R o y + 0(e 2 ). 



Thus, the function B in the denominator of z in Eq. (|29|) is given by 



B = R 2 -Q = 6RRi- l2(Rif = e 



D(3 + uj) 



mz(2t 4 ){l + u) 



+ 0(e 2 ), 



(70) 

(71) 
(72) 

(73) 



which evidently proves that B/R 2 = 0(D / z) is a small quantity justifying a posteriori 
our series expansions. In addition, since B is relatively small but finite, the term z in the 



action (29) does not diverge. 



Finally, let us determine the coefficient e a in the scale factor correction term by com- 
paring R in (71) with Eq. (40), yielding e a = — g| for dust [oj = 0) and e a = — J^ for 
radiation (w = i), respectively. Note that the term e|^ in R does not come from the 



first-order corrections oc e\2a in Eqs. (55) and (56), but it emerges as the leading term 



'4 mz 



of Ri and it provides the correct (non-zero) value for B. With these results one can check 
explicitely that z = /3(R^ 2 /^/B) m and the validity of the EOMs in Eqs. (0) and S. 



4.2 Late universe 

Analogously to the discussion on the early universe in the previous section, we determine 
now the scale factor evolution at late times, when r = p + p <C \p e \ and p e < 0. After 
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plugging B and B/R 2 from Eqs. (44) and (45) into the expressions (53) and (54) for the 
auxiliary metric h a b, we find 



1 r 



n = c 2 p e 2 1 + eei-— + 0{e z ) 

V 7 p e 

/ 4 \ — - / 3 f 
/loo = c 2 ( --mz/ej p e 2 ( 1 -ee 2 h C(e 



(74) 
(75) 



Also here, the minus sign in p e should be absorbed in the constant c 2 = \/±^-(gmz) 3 (— k) 2 , 



which drops out when calculating the connection T(h) in ( 12 ). Like before e = 1 represents 
powers of suppressed terms like B/R 2 = 0(p e /z) or r/p e , and ei j2 = 1 signal first-order 
corrections. The non-vanishing components of the connection are given by (with i = 1,2, 3) 



4 = -ee*J- + <D(e 2 ) 



r° = 


= e^ p e - + 0(e 2 , 

4mz a 


■pi 
MO " 


a 1 r ml 
~- - + eei— — + 0(e 
a 14 p e 



(76) 

(77) 
(78) 



Furthermore, we assume a constant EOS ui for the matter content as in Eq. (60), which 



simplifies the expressions for the Ricci tensor components Rq = Rq and Ri = R^ from 
Eq. (Il, 



Ro 



Ri 



a y , , r I , , . . / a 

3- + e — (1 + w)— (ei(2 + 3u) - 3e 2 - 

a 14 p e \ \a 



4mz I la/ a / 



ei- +0(e 2 ), 
a / 



(79) 



As before, the scale factor behaviour can be found by comparing the zero-order term in 
the Ricci scalar R = Rq + 3R { = 3^ + 0(e) with Eq. (4l| in the limit r -*• 0, 



k3- = p e . 

a 



(80) 



The solution is a linear combination of exponential functions, a(t) oc exp(±H e t) with 
constant .£f e . Since the shrinking solution quickly decays at late times we drop it and 
consider only the de Sitter-like solution 



(81) 



a(t) oc exp(H e t) I 1 + ee a 



where a first-order correction proportional to the constant e a has been added. Con- 
sequently, 



H = H e (l + ee a 



- 3(1+ " )r + o (e ')V ° = g ;fi +ee /' 1+3 "" 1+ "' r - =■ 



Rq = 3H 2 [l + e 



+ OW), 



3(3w - 2 + 14e„(l + 3 o;))(l + fa;)r 2 A 9ff e 2 p e 2 
h C(e ) I , Ri = e — h 0(e ), 



Upe 



Amz 



(82) 
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where e\ 2 = 1 was applied. Neglecting the tiny 0(p e /z)-teim in Ri, the Ricci scalar R 



Rq + 3Ri will be consistent with Eq. (41 ), kR = p, 



if 



p e = 3nH e and e a 

From the first condition the parameter f3 in L3 

and thus — p e 



-co 



14(l+a;) 



(83) 



4.1 



p e < can be determined because — k ~ 
—kH 2 must be of the order of the late- 
time critical energy density. This is just the effective CC corresponding to the de Sitter 
solution (81) with the tiny observed Hubble rate H e . For dust matter (uj = 0) the scale 



1/Gjv was suggested in Sec. 



factor correction oc e a vanishes and the expansion is purely de Sitter-like. 

As in the previous section, the squared Ricci tensor reads Q = (Rq) 2 + 0(e 2 ) and we 



find that the denominator of z in Eq. (29) is highly suppressed but non-zero, 



B = R 2 -Q = 6R-Ri- 12(Ri) 



, 8Egp 

2mz 



+ 0(e 2 ). 



(84) 



Therefore, our expansion in B/R 2 = 0(p e /z) turns out to be justified also at late times. 

To finalise this section we estimate the magnitude of the parameter (3 expressed as a 
power of an energy scale Mg. The dimensional analysis of z in Eq. (29) implies |/3| = Aff 
with the exponent given by d = 4 + §m. From Eq. (32) we know that z is of the order of 
the large cosmological term A, and in Eq. (83) we related (—p e ) ~ (10~ 12 GeV) 4 to the tiny 
observed energy density of the effective late-time CC. Finally, in Sec. |4.1| the parameter 
(— k) ~ (10 18 GeV) 2 was fixed by the inverse Newton constant, and now we are prepared 
to estimate f3 for given values of A and m by using p\ 



L 3 inEq. (39), 



P 



9 

2m 



„(l-m) 



(85) 



First, we note that with large values of m the tiny first factor in /3 produces small values 
of Mp = I/SI 1 ' , which makes this parameter range less attractive. On the other hand, m < 
1 provides an interesting range of energy scales. For instance, for the vacuum energy 
density |A| ~ \z\ ~ (10 16 GeV) 4 of a typical grand unified theory we obtain the following 
magnitudes of M p in units of GeV: 10~ 51 for m = 3, 1CT 2 (m = |), 10 4 (m = |), 10 12 
(m = jj), and Mg approaches z in the limit m — > 0. Apart from that, the structure of 
Eq. ( |85[ ) does not require (3 to be fixed very precisely, which constitutes a much better 
situation compared to the counterterm method in Eq. (JTl). 



5 Kottler solution 

In the previous section we have seen that the universe approaches a de Sitter cosmos 
in the limit of vanishing matter, where r — > and kR — > p e in Eq. (41). We use this 



result now for discussing the Kottler (Schwarzschild-de Sitter) solution, which describes a 
Schwarzschild black hole in the presence of a positive CC. For analysing the CC filter in 
this environment we need the 4- velocity vector u m from Eq. (26), which we know already in 



the cosmological setup. Therefore, it is useful to apply a transformation from the spatially 
flat FLRW coordinates (t, p, 9, (f>) with line element 



ds 2 



-dt 2 + a 2 {t) (dp 2 + p 2 d9 2 + p 2 sin 2 <?# 2 ) 



(86) 
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into Lemaitre coordinates (r,a,9,<f)), where the angles 6 and (f) remain untouched, 

ds 2 = -A dr 2 + A~ l da 2 + a 2 dQ 2 + a 2 sin 2 8d(p 2 . (87) 

In the latter metric the mass of the black hole in terms of the Schwarzschild radius r s and 
respectively the de Sitter radius r e define the term 



A{a) 



1 



a 



a 






The radial coordinates are related by a = a(t)p and if we require that A = 1 — (Ha) 2 with 
the Hubble rate H = -, the transformation rules between both metrics in (86) and (87) 
read 



da 
dt 



Ha = Vl- A, 



8a 

Op 



dr 
dt 



1 dr 
A' ~dp 



aHa 
~A~ 



iVT^A 
A ' 



(89) 



Next, we consider the 4- velocity vector u m in (26), which has only one non- vanishing com- 
ponent Ut = 1 in FLRW coordinates. Via the relations (89) the corresponding components 
in Lemaitre coordinates can be obtained easily, 



u T 



1, u„ 



Ha 
^4~ 



yr^A 

A ' 



(90) 



Of course, the norm u m u m 



and (90) the auxiliary metric in (26) is given by 

= ^ (9r 



h T 
K 

h T 
ii in 



1 remains invariant under this change. With Eqs. (87) 

n(-A-L u ) 



n(j, 



tt(g T 



L u u T u T ) 

— L u U T U(j) 



(91) 



nU- l -L u 



1- A 



A 2 

y/T^A 

A 



= Qgmn for the other components. 
Interestingly, h mn has non-zero off-diagonal elements, whereas the physical metric g n 



does not. Here, we have introduced the variable L Ul which follows from Eqs. (48), (52) 
and (45), and its leading term reads 

L 2 



L u '■- 



R , „ 

6-g(l + e) 



6^ 1 + e . 



(92) 



Li - L 2 a p e 

Since we consider the limit r — > 0, both terms f2 and L u are constant. Thus $7 will drop 
out from the connection T^ c according to Eq. (12). For the non- vanishing components we 
find 



Li 



TCT 



r CT 



VT^A 

1 + L u 2a A 

A + L U {A-1) 



(l + L u )2aA \a 

A + L u 



- + 2- 
a 

r s a 2 



A 2 (l-A)Tl 



-r CT 



(93) 



> L ee 



L u gyjl - A 

' {l + L u )A 



sin 



TL 



(1 + L u )2a \~a~ + T\ 



J x era 



1 

a 



1 ffd>i L <? 



cos v sm 



0, r 



l toi L ee 



cot ( 



°{A + L U ) _ 2 
' l + L u " Sm ^ 
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Moreover, the non-zero components of the Ricci tensor read 



R T 



Rn 



3(A + L U ) 
rl(l+L u y K ™~ 
3{L u (l -A) -A) 
r 2 (l + L u )A 2 ' 



Ree 



3l u vt^a 

r 2 (l + L u )A' 
3a 2 _ 
r 2 e (l+L u ) ~ 



(94) 



sin" 2 6 R r , 



Also here we find non-diagonal entries, however, they belong to the Palatini Ricci tensor, 
which is derived from h mn and not from the physical metric. Finally, we show the results 
for the scalar invariants 



R 

Q 



Rab9 



ab 



RabRc d g ac g bd 



3_ A + L u 
(R 



(95) 



2 A 



A2 +i(RoiY_ A 



A\R 



B 



R 



4 + 2L U + L 2 
(l + L u ) 2 ■ 



Q= -o 



6(2 + L U ) 
(1 + A,) 2 



R 2 



Li,, 



<7<j) 


+ 


(Ree) 2 
a 4 


+ 


sin 4 9 a 4 


(96) 


e)). 










(97) 



At leading order the last equation is consistent with Eq. (45), and it implies that B/R 2 = 
0(p e /z) <C 1 is a suitable expansion parameter. Moreover, by comparing R in Eq. ([95 



with the late-time results from Eqs. (41) and (83), we find at leading order 



R 



K 



3H 2 . 



(98) 



This implies that the de Sitter radius r e as a parameter in the Kottler metric (87) is given 
by the inverse of the final Hubble rate H e = r" 1 , just as in general relativity. Remember 
that H e originates from the small effective vacuum energy density of the order \p e \ and 
not from the large CC A. Note also that Eqs. (90) and (98) are sufficient to show that 



the metric (87) with A given in (88) is a solution of our Palatini model. The coordinate 
transformation (89) just served us to obtain the vector u m in Lemaitre coordinates. As a 



result of this section, the CC is relaxed also in situations, which can be well described by 
the Kottler metric. This can be very useful for solar system tests and the construction of 
vacuole solutions of the Einstein-Straufi type |61| [62], which we would like to discuss in 
the future. For the vacuole solutions we cannot assume r <C \p e \ as we did in Sec. 4.2, but 
the matter density r must be treated on equal footing with the effective vacuum energy 
density \p e \. 



6 Conclusions and outlook 

In the context of the old CC problem, we have presented a modified gravity model, which 
filters out a large CC A independent of its origin. Thanks to the Palatini formalism, we 
avoid problems coming from extra degrees of freedom, which are often present in the metric 
formalism. In this work, several aspects of our filter scenario have been analysed with the 
result that the standard Big Bang history of the universe is not prevented by a large CC 
term even when it dominates in size over matter or radiation. 
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Unlike many models for cosmological late-time acceleration, our setup in Eq. (29) does 
not represent a small correction to the Einstein-Hilbert term because the term z contain- 
ing the Ricci scalar R and the squared Ricci tensor Q, plays a crucial role for the filter 
effect. Thus, it is neither useful nor necessary to consider the limit z — > 0. Further- 
more, we have shown that the Hubble expansion rate in the early universe is dominated 
by the (non-vacuum) matter energy density, whereas CC contributions with equation of 
state uj = — 1 do not contribute at all in leading order. This filter effect implies a cos- 
mological background evolution similar to general relativity, where the large CC has been 
removed somehow. We have found reasonable results in the matter and radiation eras. 
In the latter the effective Newton constant is slightly increased but consistent with recent 
bounds [57J. At late times, when matter has diluted away, we enter a de Sitter phase, 
where the expansion rate H e is not dominated by the large CC, but instead it depends 



only on the magnitudes of the parameter /3 and A in Eq. (39 ). Hence, the inclusion of a CC 
counterterm which cancels A with extremely high precision is not necessary for describing 
the currently observed accelerated expansion. Our cosmological results are supported by 
the existence of a black hole solution of the Kottler type, in which the effective CC para- 
meter complies with the value of the final de Sitter expansion rate H e . This suggests that 
cosmology and the astro-physical domain can be smoothly connected. 

The robustness of the CC filter can be seen also from a different perspective. So far we 
have considered A to be a large constant generally of the order of the largest contribution 
to the CC. During the cosmic evolution, however, it is not unlikely that shifts <5A of 
vacuum energy occur, as a result of phase transitions for instance. We have shown that 
these contributions, which are smaller than A by definition, have no influence on gravity at 



leading order if we treat them as part of the matter sector described by r = p+p in Eq. ( 37 ) . 
Alternatively, one could shift the cosmological term by A — > A + 5A, but also this procedure 
does not change our results because the differences are of the order |<5A/A| <C 1 contributing 
only to the small correction terms. Hence, the exact value of A is not important for the CC 
filter effect even when vacuum shifts 5 A arise dynamically as the universe evolves. Note 
that 5A can be much larger than the current critical energy density despite being small 
compared to A. 

Let us now sketch some open points and discuss possible generalisations of our setup, 
which we want to address in the future. First, it is important to mention that our filter 
effect is based on the largeness of the CC and its universal energy-momentum structure. 



This can be seen nicely in Eq. (34), where the CC contribution to the Ricci scalar R can 



be completely removed by a suitable model building construction, which leads to Eq. (37). 



It is clear that as long as A dominates over other sources in the energy-momentum tensor 



and its trace (31 ), the procedure will be the same, which removes from R the dominating 



term z ~ A as given in Eq. (32 ). Therefore, we can at least conjecture that matter sources 
different from the perfect fluid form as well as other background metrics exhibit also the 
CC filter property. This would be useful for future studies of the Newtonian limit as well 
as perturbations in cosmology and astro-physical setups, which are more involved than the 
Kottler solution. For the latter, it seems that vacuole space-times can be easily constructed 
once a solution is found which interpolates between the matter era and the final de Sitter 
phase. The smoothness of R as shown in Fig. [T] suggests such a solution. More difficult 
could be the analysis of the very early cosmological epoch when matter dominated over A 
(if there was such a time) or during primordial inflation, respectively. In that case many 
approximations we used for the subsequent eras cannot be applied anymore, and a new 
analysis is necessary. 



18 



Apart from that, we have shown only one numerical example to support our analytical 
considerations, and it would be nice to have more numerical results discussing e.g. the 
transition phases between different epochs and the constraints on the parameters (3 and m. 



Also the L3 > case in Eq. (39) might correspond to an interesting solution, maybe it 



describes a contracting universe at late times because a vanishing matter density is not a 



solution of Eq. (39). However, since we observe an accelerating cosmos, the L3 < case 
discussed here is preferred. 

Moreover, we should mention that Palatini models are subject to constraints from 
astrophysical bodies and even small scale atomic physics, see e.g. |63|, [35| WI\ l4"8l 137] . 
However, many results correspond to /(i?)-type actions or to models, where the Einstein- 
Hilbert term is amended by a correction, which becomes large for low curvature. It is clear 
that our setup is not of this type. We will see in the future whether this is an advantage 



or not. Nevertheless, it might help that the Ricci scalar in Eq. (41) and the denominator 



function B in (44) approach a non-zero constant in the limit of vanishing matter. This 
suggests a stable vacuum, which limits the values of the geometrical scalars from below. 
Finally, it has been argued in Refs. |47| HH] that the averaging over matter sources in the 
Palatini framework has some non-trivial aspects, too. 

In the end, it would be interesting to know whether the CC filter effect is just a curiosity 



of our specific model in Eq. (29), or if there are completely different choices, which have 



the same property. We want address these questions in the future. 
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